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Effect of shear viscosity on elliptic flow is studied in causal dissipative hydrodynamics in 2+1 
dimensions. Elliptic flow is reduced in viscous dynamics. Causal evolution of minimally viscous 
fluid (77/s=0.08), can explain the PHENIX data on elliptic flow in 16-23% Au+Au collisions up to 
Pt ~3.6 GeV. In contrast, ideal hydrodynamics, can explain the same data only up to pr «1.5 GeV. 
Pt spectra of identified particles are also better explained in minimally viscous fluid than in ideal 
dynamics. However, saturation of elliptic flow at large pr is not reproduced. 
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One of the important results in Au+Au collisions is 
the large elliptic flow is non-central collisions [j], 0, H, 0] • 
Large elliptic flows establish that in non-central Au+Au 
collisions, a collective QCD matter is created. Whether 
the matter can be characterized as the lattice QCD Q 
predicted Quark-Gluon-Plasma (QGP) or not, is still a 
question of debate. Qualitatively, elliptic flow are nat- 
urally explained in hydrodynamics, rescattering of sec- 
ondaries generates pressure and drives the subsequent 
collective motion. In non-central collisions, the reaction 
zone is asymmetric (almond shaped), pressure gradient 
is large in one direction and small in the other. The 
asymmetric pressure gradient generates the elliptic flow. 
As the fluid evolve and expands, asymmetry in the re- 
action zone decreases and comes a stage when reaction 
zone become symmetric and system no longer generate 
elliptic flow. Elliptic flow is early time phenomena and a 
sensitive probe to the early stage of the fluid. 

Ideal hydrodynamics has been quite successful in ex- 
plaining RHIC data on elliptic flow [6|. Assuming that 
the collision produces QGP, which undergoes 1st order 
phase transition at a critical temperature T c =164 MeV, 
elliptic flow of identified particles are explained up to 
Pt ~1.5 GeV. Ideal hydrodynamics also reproduces the 
transverse momentum dependence of identified particles, 
again up to pt ~ 1.5 GeV, which exhaust more than 
99% of particle production. Success of ideal hydrody- 
namics in explaining bulk of the data Q , together with 
the string theory motivated lower limit of shear viscosity 
v/ s > 1/47T @, H[ has led to a paradigm [1, [13] that in 
Au+Au collisions, a nearly perfect fluid is created. 

However, the paradigm of perfect fluid, produced in 
Au+Au collisions at RHIC, need to be clarified. As in- 
dicated above, the ideal hydrodynamics is only part ially 
successful and in a limited px range (jjt <L5 GeV) [ill ]. 
The transverse momentum spectra of identified particles 
starts to deviate form ideal fluid dynamics prediction be- 
yond pt ~ 1.5 GeV. Experimentally determined HBT 
radii are not reproduced in the ideal fluid dynamic mod- 
els, the famous "HBT puzzle" [12]. It also does not repro- 
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duce the experimental trend that elliptic flow saturates 
at large transverse momentum. Contrary to experiment, 
ideal hydrodynamics predicts continually increasing el- 
liptic flow. These shortcomings of ideal fluid dynamics 
indicate greater importance of non-equilibrium (dissipa- 
tive) effects in the pt ranges greater than 1.5 GeV. In- 
deed, saturation of elliptic flow at large is a manifesta- 
tion of non-equilibrium effects. In parton cascade simu- 
lations of relativistic Boltzmann equation, where dissipa- 
tive effects are naturally included, elliptic flow saturates 
at large pt [131 ]. However, parton cascade simulations do 
not reproduce experimental data unless the elastic par- 
ton cross-section is very high, a ~ 45 mb. Realistic cross- 
sections do not generate enough flow. Recently, elliptic 
flow is studied in a hybrid " hydro+cascade" model [l4| . 
In the model, initially produced QGP evolve following 
the ideal hydrodynamics. Just below the phase transi- 
tion, a hadronic transport model is employed for the late 
stage evolution. The model seems to require additional 
dissipation at the early quark-gluon plasma phase. 

Dissipation in the early QGP phase can be conve- 
niently studied in relativistic dissipative hydrodynam- 
ics. Theories of relativistic dissipative hydrodynamics 
are well developed [HI [H, [llj • Problem of causality 
violation in 1st order theories [HI [l6| . are removed in 
the Israel-Stewart's [TtJ causal, 2nd order theories. Re- 
cently, there has been significant progress in numeri- 
cal implementation relativistic dissipative hydrodynam- 
ics liliiiiiaiMISIilliilSlSl^lS- Ellip- 
tic flow in causal hydrodynamics has been calculated in 
[27L [28| . Romatschke and Romatschke 127] concluded 
that while data on integrated «2 is consistent with a ra- 
tio of viscosity over entropy density up to r)/s «0.16, the 
data on minimum bias vi favor a much smaller viscosity, 
rj/s «0.03, lower than the ADS/CFT bound ?//s=0.08. A 
different result was obtained in [28|. They did not com- 
pare with experiments but minimal viscosity 77/5=0. 08 
leads to large reduction in elliptic flow compared to 
ideal hydrodynamics. At the Cyclotron Centre, Kolkata, 
we have developed a code, " AZHYDRO-KOLKATA" to 
solve '"causal"' dissipative hydrodynamics. In an ear- 
lier publication [26|, we have discussed the space-time 
evolution of QGP fluid without any phase transition. 
AZHYDRO-KOLKATA results for elliptic flow (includ- 
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ing phase transition) are presented here. 

Elliptic flow of a particle (say ir~) is quantified as the 
2nd harmonic of azimuthal distribution of 7r~, 



V2(pt) 



r d 2 N 

dyd' 2 pT 



(i) 



where dyd^p T ^ s ^ ne mvar i an t distribution of n . In 
Cooper-Frye prescription, the invariant distribution is 
obtained as, 



dN 
dyd 2 p q 



dZ^f(x,p), 



(2) 



where S M is the freeze-out hyper surface and f(x,p) is the 
distribution function. In viscous dynamics, the system is 
not in (local) equilibrium. If the system is not far from 
equilibrium, f(x,p) can be approximated as, 



Eqg]). 



on the other hand grow less than linearly 



with px- An impportant conclusion can be reached from 
EqEl asympotically, viscous hydrodynamics do not pre- 
dict saturation of elliptic flow. 
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Au+Au@b=0 at t=8 fm 

dx=0.2,dx=0.02 

dx=0.2,dx=0.01 

dx=0.1,dx=0.01 
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f(x,p) = f^(x,p)[l + cf>(x,p)], 



(3) 



where f°(x,p) is the equilibrium distribution function, 
and <p(x,p) << 1, is a small correction due to non- 
equilibrium effects. With shear viscosity as the only dis- 
sipative force the correction factor (j)(x,p) can be locally 
approximated as, 



(j){x,p) 



e^p v 



2T 2 (e+p) 



Wp". (4) 



where e,p and T axe the (local) energy density, pressure 
and temperature. 7r MI/ 's are the shear stress tensors. 

Accordingly, in viscous dynamics, n~ invariant distri- 
bution has two parts, 



dN 



dyd 2 pT 



d^f°(x,p) + J dZ^[f(x,p)ct>(x,p)} 
dN e i dN ne i 

(5) 



dyd 2 pT dyd 2 pT 



where d y d z pT is the non-equilibrium correction to the w 
equilibrium distribution d y d 2 pT ■ Since (f){x,p) << 1 it is 



then necessary that dyd 2 pi 
to the 1st order, 
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Expanding Eq[T] 



V2{pt) = vl q {pT) 
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where v 2 is the equilibrium elliptic flow and V2° rr is 
the correction due to the non-equilibrium effects. Non- 
equilibrium correction grow quadratically with pt (see 



FIG. 1: (color online) constant energy density contours in 
b=0 Au+Au collisions after 8 fm of evolution. 

Assuming longitudinal boost-invariance, we have 
solved causal dissipative hydrodynamics in 2+1 dimen- 
sion. We have included the shear viscosity only. Details 
of equation solved and some results for space-time evo- 
lution could be found in (25[. To simulate the Au+Au 
collisions at RHIC, we have used the standard initial con- 
ditions described in Q. At initial time Tj=0.6 fm, the 
QGP fluid was initialized with central entropy density 
Sini = 110/to~ 3 . It corresponds to peak energy density 
~ 35GeV/fm 3 . We also use the equation of state EOS-Q 
described in Q , incorporating a first order phase transi- 
tion at a critical temperature T c =164 MeV. In viscous hy- 
drodynamics, additionally, one has to initialize the shear 
stress-tensors. We assume that at initial time r^, shear 
stress-tensors have attained their boost-invariant values 
[25| . Viscous hydrodynamics introduces some new pa- 
rameters, the viscosity coefficient 77 and the relaxation 
time T re i x . Viscosity coefficient of a QGP or hadronic 
fluid is quite uncertain. In a strongly coupled field the- 
ory, ADS/CFT estimate gives a lower limit to shear vis- 
cosity, 2 > _L ps 0.08 @, Presently we show results 
for the minimal shear viscosity, ?7/s=0.08. For the relax- 
ation time, we use the kinetic theory approximation for 
a Boltzmann gas, T reix = ^ « 

It is very important to establish the accuracy of the 
numerical code. Being a ratio, elliptic flow is a sensitive 
observable and any small error in computation can lead 
to large uncertainty in the elliptic flow. Viscous hydro- 
dynamics do not posses analytical solution against which 
the numerical results could be tested. We use the general 
procedures to test a numerical code, (i) the results should 
be stable against change in integration step lengths, (ii) 
any symmetry in the system should not be destroyed and 
(iii) unphysical maxima or minima should not creep into 
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due to numerical error. In Fig[TJ we have shown the 
constant energy density contours in a zero impact pa- 
rameter Au+Au collisions after 8 fm of evolution. The 
contours are drawn for three different sets of integration 
step lengths, (ijdx = dy = 0.2 fm, dr = 0.02/m, (the 
black lines) , (ii)dx = dy = 0.2 fm, dr = 0.01/m, (the 
red lines) and (m)dx = dy — 0.1/m, dr — 0.01/m, (the 
blue lines). The ratio of viscosity over entropy density is 
?7/s=0.08. FigfT]shows that all the three conditions men- 
tioned above are satisfied. The evolution is stable against 
change in integration step lengths. In the interior of the 
fluid, energy density surfaces are exactly reproduced and 
one can not distinguish the lines for different integration 
step lengths. At large radius, energy density distribution 
does show dependence on integration step lengths but 
the dependence is small. For example, the spatial uncer- 
tainty in the e=0.1 GeV surface is less than 3-4% only. 
In a b=0 collision, initial energy density is azimuthally 
symmetric. The symmetry is maintained during the evo- 
lution. Also the smooth fall of energy density along the 
radius indicates that no new maxima or minima are gen- 
erated during the evolution. 
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FIG. 2: (color online) The solid and dashed lines show the pr 
dependence of elliptic flow in evolution with integration step 
lengths dx = dy — 0.2, dr = 0.02 and dx — dy = 0.1, dr = 
0.01. The black lines corresponds to ideal fluid. The red and 
blue lines corresponds to viscous fluid with ?7/s=0.04 and with 
n/s=0.08 respectively. The filled circles are PHENIX data 
[29l | on elliptic flow of charged particles in 16-23% Au+Au 
collisions. 

Any dependence of the hydrodynamic evolution on in- 
tegration step lengths will be reflected on the freeze-out 
surface and consequently on observables like particle px 
spectra, elliptic flow etc. In FigfJl we have computed 
the elliptic flow for ir~ in a b=6.5 fm Au+Au collision. 
The dashed and solid lines are V2 in evolution with inte- 
gration step lengths (i) dx = dy = 0.1/m, dr — 0.01/m 
and (ii) dx = dy — 0.2/m, dr = 0.02 fm respectively. 
We have shown results both for ideal (black lines) and 
viscous fluid with two values of viscosity ?y/s=0.04 (the 
red lines) and 0.08 (the blue lines). The freeze-out tem- 



perature is assumed to be Tp=150 MeV. Both in ideal 
and in viscous dynamics, i>2 depend marginally on in- 
tegration step lenghs. Elliptic flow in evolution with 
dx = dy — 0.1/m, dr — 0.01/m is ~ 10% more than in 
evolution with dx — dy — 0.2 fm, dr = 0.02/m. Elliptic 
flow is very sensitive observable. Stablity of V2 imply that 
the hydrodynamic evolution, whether computed with in- 
tegration step lengths (i) dx = dy = 0.1/m, dr = 0.01 fm 
or (ii) dx = dy = 0.2/m, dr — 0.02 fm, lead to nearly 
identical freeze-out surfaces. 
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FIG. 3: pr dependence of elliptic flow in ideal (the black 
lines) and in minimal viscous (the blue lines) dynamics for 
T F =130,140,150 and 160 MeV in a b=6.5 fm Au+Au collision. 
The filled circles are the PHENIX data [1| on the elliptic flow 
in 16-23% Au+Au collisions. 
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FIG. 4: PHENIX data [33] on the transverse momentum de- 
pendence tv~ , K + and proton in 0-5% Au+Au collisions are 
shown. The blue and black lines are predictions from mini- 
mally viscous hydrodynamics and from ideal hydrodynamics 
in a b=2.3 fm Au+Au collision. 

Effect of viscosity on elliptic flow is also evident in 
FigEl Under identical conditions elliptic flow is less 
in viscous dynamics than in ideal dynamics. More vis- 
cous is the fluid less is the elliptic flow. b=6.5 fm col- 
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lision roughly corresponds to 13-26% Au+Au collision. 
In FigH PHENIX data [H[ on the elliptic flow in 16- 
23% Au+Au collisions are shown (the filled circles). As 
mentioned earlier, ideal fluid dynamics fails to explain 
the data beyond pt «1.5 GeV. Data are not explained 
in minimally viscous dynamics also. ADS/CFT lower 
bound on viscosity 7y/s=0.08 reduce the elliptic flow more 
than necessary. Reduction is less for still lower viscosity 
and PHENIX data are well explained up to pr ~3 GeV 
for viscosity ?7/s=0.04. A similar result is also obtained 
in ■ Minimum bias STAR data seems to favor smaller 
viscosity (?7/s=0.03), lower than the ADS/CFT bound. 

In a viscous fluid, elliptic flow depend strongly on 
the freeze-out temperature. For the minimally viscous 
fluid (r?/s=0.08), PHENIX data could be explained if 
the hadronic freeze-out occur at temperature lower than 
150 MeV. The shear stress tensors will evolve for longer 
duration to lower values and contribute less to the non- 
equilibrium distribution function and to the elliptic flow. 
In FiglJl the elliptic flow for freeze-out temperatures 
7>=130,140,150 and 160 MeV are shown. Curves are 
drawn up to pr such that 

< jyv^i As expected, 
elliptic flow increases as Tp decreases, pr range over 
which viscous dynamics remain valid is also extended. 
For Tp=130 MeV, minimally viscous dynamics reproduce 
the PHENIX data in 16-23% Au+Au collisions, well up 
topT=3.6 GeV. However, viscous dynamics do not repro- 
duce saturaton. At large pr jets are important. With- 
out the effect of jets accounted for, possibly saturation 
of elliptic flow could not be explained. For comparison, 



we have also shown the elliptic flow in ideal dynamics. 
In ideal dynamics freeze-out dependence is compartively 
less. We may mention that fit to elliptic flow data in other 
centrality ranges of collisions is not as good as shown in 
FigH Minimally viscous dynamics produces less elliptic 
flow in very central collisions (0-5%, 5-10%) collisions. 
In less central collisions data are reasonably well repro- 
duced. 

Viscous hydrodynamics with minimal viscosity 
77/s=0.08 also reproduces the transverse momentum 
distribution of identified par ticles. In FigfJ] we have 
shown the PHENIX data [30| on transverse momentum 
distribution of 7r~, K + and protons in 0-5% Au+Au 
collisions. The blue lines are viscous hydrodynamics 
predictions in a b=2.3 fm Au+Au collision. Freeze- 
outtemperature is 7V=130 MeV. Minimally viscous 
hydrodynamics (?7/s=0.08) reproduces the 7r~, K + and 
proton spectra up to j?t=3 GeV. A comparable fit could 
not be obtained in ideal dynamics. 

To summarise, in causal dissipative hydrodynamics, we 
have studied the effect of viscosity on elliptic flow. Un- 
der identical conditions, elliptic flow is less in viscous 
dynamics than in ideal dynamics. For minimally viscous 
fluid (?7/s=0.08), PHENIX data on elliptic flow in 16- 
23% Au+Au collisions are well explained up to pr «3.6 
GeV. This is to be contrasted with ideal dynamics which 
can explain the flow up to px ~ 1.5 GeV only. Minimally 
viscous fluid also explains the px spectra of 7r~ , K + and 
protons, much better than in ideal dynamics. 
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